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ABSTRACT 



Context. The observed abundance peculiarities of many chemical species relative to the expected cluster metallicity in 
blue horizontal-branch (BHB) stars presumably appear as a result of atomic diffusion in the photosphere. The slow 
■ rotation (typically vsini < 10 km s~^) of BHB stars with effective temperatures Tcs > 11,500 K supports this idea 

$ZL(' since the diffusion mechanism is only effective in a stable stellar atmosphere. 

Aims. In this work we search for observational evidence of vertical chemical stratification in the atmospheres of six hot 
BHB stars: B84, B267 and B279 in M15 and WF2-2541, WF4-3085 and WF4-3485 in M13. 

Methods. We undertake an abundance stratification analysis of the stellar atmospheres of the aforementioned stars. 



o 



o 



- 

X 



^ ' based on acquired Keck HIRES spectra. 

' Results. We have found from our numerical simulations that three stars (B267, B279 and WF2-2541) show clear 
signatures of the vertical stratification of iron whose abundance increases toward the lower atmosphere, while the 

^SJ ■ other two stars (B84 and WF4-3485) do not. For WF4-3085 the iron stratification results are inconclusive. B267 also 

^ ■ shows a signature of titanium stratification. Our estimates for radial velocity, vsini and overall iron, titanium and 

' phosphorus abundances agree with previously published data for these stars after taking the measurement errors into 

I account. The results support the hypothesis regarding the efficiency of atomic diffusion in the stellar atmospheres of 

! BHB stars with T^s > 11,500 K. 

Key words, stars: atmospheres - stars: horizontal-branch - stars: chemically peculiar 

1. Introduction The Hertzsprung-Russell diagrams of some globular 

clusters show long blue tails (an extension of the HB), 

According to the current understanding of stellar evolution, populated by very hot BHB stars and extreme horizontal 

the horizontal-branch (HB) stars are post-mam sequence ^^.^nch (EHB) stars. PubUshcd data on the BHB stars ar- 

^ stars that burn hchum m their core and hydrogen m a shell ^hat the hot BHB stars show remarkable difTerences in 

H , (e.g. Moehler,2004)- In this paper we consider the HB stars physical properties when compared to the cool BHB stars. 

. - . that are located m the blue part of the HB, to the left of Using high-precision photometry of stars in M13, Ferraro 

the RR Lyrae mstability strip. Most researchers call them ^1. (HMHD have found gaps in the distribution of stars 

blue horizontal-branch (BHB) stars to distinguish them ^long the blue tail. One of these gaps, labeled as Gl, is 

from the red horizontal-branch (RHB) stars, which exhibit located at T,s -11,000-12,000 K. Grundahl et al. (|T998l) 

different observational properties. Sandage & Wallerstem ^^^^ ^^e results of Stromgren ww6y/3-photometry finding 

UMH have found from analysis of the colour-magnitude g^od agreement between the theoretical prediction of stel- 

diagrams of globular cluster^ that the HB generally be- ^^r evolution models and the observed location of BHB 

comes bluer with decreasing metallicity Derived masses of g^ars, except for the hot BHB stars, whose w-magnitudes 

the cool (Teff <11,500 K) BHB s tars in the globular cluster ^re brighter than predicted. It appears that this u-jump 

NGC 6388 (Moehler & SweigartEDDB]) are m a good agree- jg observed for the hot BHB stars and coincides with the 

ment with the predictions of canonical HB evolution, except temperature range of the Gl gap in M13. Similar u-jumps 

for the hot BHB stars with T^s >11,500 K, where the es- ^ave also been found for other globular clusters (Grundahl 

timated stellar masses seem to be lower than the canonical ^1. [T9991) . For the hot BHB with effective temperatures 

values. yp ^o 20,000 K, the surface gravities derived from the fits 

■ of Balmer and helium line profiles appear to be lower than 

Send offprint requests to: V. Khalack predictions of stellar evolution models, while the gravi- 
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, 1 • i T^r-i noofir, lu good agreement With thcoretical prcdictious (Mochler ct 

Avenue NW, Washington DC 20392 al. IMS. iMIa, IMZfl IlQOl. The stellar rotation velocity 

Most of the known BHB stars are found m globular clusters. * 
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Table 1. Journal of Keck+HIRES spectroscopic observations of the selected hot BHB stars from Behr (j2003b[) . 



Cluster/Star 


HJD 


Exposure 


S/N 


Coverage 


Seeing 


Slit' width 




2450000+ 


Time (s) 




(A) 


(arcsec) 


(arcsec) 


M13/WF2-2541 


1046.7902 


3x1500 


44 


3885-6292 


0.90 


0.86 


M13/WF4-3085 


1052.7338 


3x1200 


37 


3888-5356 


1.10 


0.86 


M13/WF4-3485 


1053.7793 


3x1200 


34 


3888-5356 


0.90 


0.86 


M15/B84 


1053.9024 


4x1400 


34 


3888-5356 


0.80 


0.86 


M15/B267 


1053.9731 


4x1400 


26 


3888-5356 


0.80 


0.86 


M15/B279 


1053.8312 


4x1400 


34 


3888-5356 


0.90 


0.86 



distribution of BHB stars also appears to have a discon- 
tin uity at T^s - 11,500 K (Pe terson et al. fTOM l Behr et 
al. I2000al Recio-Blanco et al. I2004p . indicating that the 
hotter stars show modest rotation with wsinz < 10 km 
s~^, while the cooler stars rotate more rapidly on average. 
Comprehensive surveys of abundances also show that the 
hot BHB stars have abundance anomalies when compared 
to the cool BHB stars in the same cluster (Glaspey et al. 
[Ml Grundahl et al . [TW^l Behr et a l. [TOMl [20n0bl Behr 
I2003al Fabbian et al. [20051 Pace et al. l2006l) . 



The observed phenomena such as the low gravity, pho- 
tometric jumps and gaps, abundance anomalies and slow 
rotation suggest that atomic diffusion could be important 
in the stellar atmospheres of hot BHB stars. Atomic diffu- 
sion arises from the competition between radiative accel- 
eration and gravitational settling. This can produce a net 
acceleration on atoms and ions, which results in their diffu- 
sion through the atmosphere (Michaud JL970) . In order for 
atomic diffusion to produce a vertical stratification of the 
abundances of particular elements, the stellar atmosphere 
must be hydrodynamically stable. According to Landstreet 
(|1998p . photospheric convection should be very weak at 
the effective temperatures of BHB stars. Theoretical at- 
mospheric models of Hui-Bon-Hoa, LeBlanc & Hauschildt 
(|2000p showed that the observed photometric jumps and 
gaps for hot BHB stars can be explained by elemental dif- 
fusion in their atmosphere. Behr (2003b]) has shown that 
adoption of a microturbulent velocity of or 1 km s~^ pro- 
vides the best fit to line strengths in the spectra of hot BHB 
stars. This fact supports the proposal that strong velocity 
fields are not present in the atmospheres of hot BHB stars. 



While synthesizing spectral line profiles, Khalack et al. 
(|2007p have recently found vertical abundance stratifica- 
tion of sulfur in the atmosphere of the field BHB star 
HD 135485. In this paper we also attempt to detect signa- 
tures of vertical abundance stratification of elements from 
line profile analyses of several other BHB stars for which we 
have appropriate data. Together with the data on stratifi- 
cation of the sulfur abundance in HD 135485, new positive 
results would provide a convincing argument in favour of 
efficient atomic diffusion in the atmospheres of hot BHB 
stars. In Sec. [5] we discuss the properties of the acquired 
spectra, while in Sec. [3] we describe details concerning the 
simulation routine and adopted atmospheric parameters for 
the program stars. The evidence for vertical stratification 
of some chemical species is given in Sec. [31 while the esti- 
mation of mean abundances and velocities is described in 
Sec. [5] A discussion follows in Sec. [51 



2. Observations 

In this paper we have selected hot BHB stars from the list 
of objects published by Behr Ij2003bp and which have com- 
paratively high signal-to-noise ratio (S/N) spectra avail- 
able. Spectroscopic observations of the selected stars were 
undertaken in August 1998 with the Keck I telescope and 
the HIRES spectrograph. The journal of spectroscopic ob- 
servations is shown in Table [1] where individual columns 
give the object identification, the heliocentric Julian Date 
of the observation, the exposure time, the S/N per pixel the 
spectral coverage, the size of the seeing disk (FWHM) and 
the size of the CI slit that provides a spectral resolution 
oi R = X/SX =45,000. For the aforementioned stars, Behr 
(j2003bp has found that the underfilling of the slit should 
not change the estimated spectral resolution by more than 
4%-7%. 

The package of routines developed by McCarthy l|1990p 
for the FIGARO data analysis package (Shortridge 1993|) 
was employed to process the spectra. A comprehensive de- 
scription of the data acquisition and reduction procedure 
is presented by Behr (|2003bp . This followed a standard 
prescription with bias subtraction, flat-fielding, order ex- 
traction, and wavelength calibration from thorium-argon 
comparison lamp observations. To minimize the potential 
distortion of narrow spectra features, cosmic ray hits were 
identified and removed by hand. 

3. Line profile simulations 

3.1. Stellar atmosphere parameters 

The line profile simulations were performed with the 
Zeeman2 spectrum synthesis code (Landstreet 119881 
Wade et al. l200ip . The stellar atmosphere models were cal- 
culated with the Phoenix code (Hauschildt et al. I1997P as- 
suming LTE (Local Thermodynamic Equilibrium) and us- 
ing the stellar atmosphere parameters extracted from Behr 
(|2003bp and listed in Table [H In calculating the stellar 
atmosphere models for program stars, we have used so- 
lar metallicity with the enhanced iron and depleted helium 
abundances derived by Behr (12003bp . The depleted helium 
abundance has also been taken into account during line pro- 
file simulations using Zeeman2. We used gaussian instru- 
mental profiles with widths derived from the comparison 
arc spectra. 

To simulate the spectra of BHB stars, Behr (j2003bp ini- 
tially adopted a microturbulent velocity ^ =2 km s~^ and 
subsequently updated it during the simulation routine ac- 
cording to the results of his model fits. He has obtained best 
fit values of ^ from to 0.9 km s^^ for the BHB stars stud- 
ied here (see Table [2]) . Since we search for possible vertical 
abundance stratification, our results can be affected by mi- 
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logTsooo logxsooo 

Fig. 1. Abundance estimates from the analysis of (a) Fell and (b) Tin lines as a function of line (core) formation 
optical depth assuming ^=0 km s^^ (filled circles) and ^=2 km s^^ (open circles) for B267. The solid and dashed lines 
approximate the filled and open circles respectively with straight line using the least-square algorithm. 



a) b) 




logTsooo logTsooo 

Fig. 2. The same as in Fig. [U but for (a) Fell and (b) Tin lines extracted from the B279 spectrum. 



Table 2. Parameters of stellar atmospheres for the selected 
hot BHB stars from Behr (|2003bp . 



Cluster/Star 


Teff 


logs 




K 




(K) 


(dex) 


(km s"^) 


(km s'^) 


M13/WF2-2541 


13000 


4.0 


0.0 


-257.5 


M13/WF4-3085 


14000 


4.0 


0.0 


-255.7 


M13/WF4-3485 


13000 


4.0 


0.0 


-246.8 


M15/B84 


12000 


3.5 


0.5 


-108.2 


M15/B267 


11000 


3.5 


0.0 


-114.4 


M15/B279 


11000 


3.5 


0.9 


-104.4 



croturbulence. To estimate this influence we have adopted 
here two different microturbulent velocities ^ =0 and 2 km 
s~^ for our simulations. 

3.2. Procedure 

We examined the spectrum of each star to establish a 
line list suitable for abundance and stratification analysis. 



The line identification was performed using the VALD-2 
(Kupka et al. 119991 Ryabchikova et al. 1999), Pickering et 
al. (|2001l) and Raassen & Uylingfl (|1998)) line databases. 
For all of the program stars the best represented element 
(with a number of readily visible line profiles) is iron. Some 
stars also show Tin or Pn lines. Therefore iron, titanium 
and phosphorus were selected for analysis of their possi- 
ble vertical abundance stratification. For Fen lines, we ex- 
tracted atomic data from Raassen & Uylings (1998), while 
for Pn lines we used VALD-2. Atomic data for Tin lines 
were taken from Pickering et al. (|200ip . 

The Zeeman2 code has been modified 
(Khalack & Wade I2006P to allow for an automatic 
minimization of the model parameters using the downhill 
simplex method (Press et al. I1992p . The minimization 
routine finds the global minimum of the function, which 
is specified as the measure of differences between simulated 
and observed line profiles. The relatively poor efficiency of 
the downhill simplex method, requiring a large number of 

^ I ftp://ftp.wins.uva.nl/pub/orthl 



4 



Khalack et al.: Search for vertical stratification of metals in atmospheres of BHB stars 
a) b) 




0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 
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Fig. 3. Dependence of iron abundance on the lower level excitation potential of Fe ii lines extracted from (a) B267 and 
(b) WF4-3485 spectra. The linear least-square fit to the data is shown by solid lines. Here the star B267 shows direct 
evidence of iron stratification, while WF4-3485 does not. 



Table 3. List of iron spectral lines used for the abundance analysis. 



A, A 


logs/ 




log 7rad 


Ref. 




Object 


name 






Fen 


WF2-2541 


WF4-3085 


WF2-3585 


B84 


B267 B279 


4122.668 


-3.300 


20830.58 


8.49 


R&U 




X 




X 




4128.748 


-3.578 


20830.58 


8.61 


R&U 




X 






X 


4173.461 


-2.617 


20830.58 


8.61 


R&U 


X 


X 


X 


X 


X 


4177.618 


-3.776 


73395.93 


8.61 


R&U 




X 




X 


X 


4177.692 


-3.449 


20516.96 


8.61 


R&U 




X 




X 


X 


4178.862 


-2.535 


20830.58 


8.49 


R&U 


X 


X 


X 


X 


X 


4184.261 


-1.938 


90397.87 


8.35 


R&U 




X 









Note 1. Table [3] is presented in full in the electronic edition of Astronomy and Astrophysics. A portion is shown here for guidance 
and content. 



function evaluations, is a well known problem. However, 
repeating the minimization routine several times in the 
vicinity of a supposed minimum in the chosen parameter 
space allows us to verify if the method converges to the 
global minimum (for more details about this minimization 
routine see Khalack et al. I2007p . 

To search for the presence of abundance stratification, 
we have estimated the abundance of a chemical element 
from an independent analysis of each selected line profile. 
This method operates with three free model parameters 
(the element's abundance, the radial velocity K- and wsini) 
that are derived from each line profile using the aforemen- 
tioned automatic minimization routine. To analyse the ver- 
tical abundance stratification in this method we build the 
scale of optical depths T5000 for the list of selected line pro- 
files. First, we calculate the line optical depth in the line 
core for every layer of the stellar atmosphere model. Next, 
we suppose that each analyzed profile is formed mainly at 
T(,=l, which corresponds to a particular layer of the stellar 
atmosphere with respective continuum optical depth T5000 ■ 
All simulations are performed with stellar atmosphere mod- 
els that contain 50 layers. 

In general, we have selected for our analysis lines found 
to be free of predicted or inferred blends. However, if a 
blend is from a line of the same element that forms the main 
line profile, such a line was also included in our simulation. 
When a simulated line profile results in a radial velocity 



which differs significantly (more than 2 km s"^) from the 
average for the analysed star, we exclude such a line 
from further consideration. The difference in radial velocity 
may be evidence of line misidentifications or inaccurate line 
wavelengths. Atomic data (and sources) for the final list of 
iron lines selected for our study is given in Table [3l 

4. Vertical abundance stratification 

Applying the technique described in Subsec. 13.21 we have 
attempted to determine if the iron abundance is vertically 
stratified in the atmospheres of the six selected BHB stars. 
For comparison, we have also investigated the possible ver- 
tical stratification of titanium and phosphorus for some 
BHB stars, where these elements are represented by a suf- 
ficient number of spectral lines. 

All of the selected BHB stars are slowly rotating and 
have sharp absorption lines that result in only 8 to 13 spec- 
tral bins per line profile. This fact, together with the S/N 
and the uncertainties in the atomic data, are the main con- 
tributors to errors in the abundance inferred from a single 
line. This makes the detection of weak vertical abundance 
variations difficult. 

Our abundance analysis for ^ = kms""'^ shows that in 
the atmospheres of four BHB stars (B267, B279, WF4-3085 
and WF2-2541) the iron abundance generally increases to- 
wards the deeper atmospheric depths (see Figs. [T] to [5]). It 
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Fig. 4. The same as in Fig. [1] but for (a) Fe ii and (b) P ii lines extracted from the WF4-3085 spectrum. Here the visual 
strong rise in iron abundance at large optical depth manifests itself by four lines and may be real. Nevertheless, to draw 
the final conclusion we need more data for this range of optical depths. 




Table 4. Slopes of vertical abundance stratification calculated for the linear least-square fit to the data assuming zero 
microturbulence. 



Cluster/Star 


Slopes of \og{Fe/Ntot) vs. 


n 


Slopes of log 


[P/Ntot) vs. 


n 




logE^ 


log rsooo 




logs. 


log rsooo 




M13/WF2-2541 


0.18±0.04 


0.17±0.03 


53 


-0.67±0.31 


-0.14±0.06 


16 


M13/WF4-3085 


0.10±0.03 


0.15±0.02 


117 


0.54±0.31 


0.24±0.04 


12 


M13/WF4-3485 


-0.06±0.03 


-0.02±0.03 


50 








M15/B84 


0.04±0.05 


0.03±0.04 


25 
















Slopes of log{Ti/Ntot) vs. 




M15/B267 


0.49±0.06 


0.33±0.03 


46 


0.47±0.10 


0.36±0.07 


16 


M15/B279 


0.23±0.04 


0.23±0.02 


71 


0.21±0.11 


0.30±0.08 


11 



should be noted that the range of optical depths diagnosed 
by the iron lines is much smaller in the stars B279 and WF2- 
2541 than in the other two stars shown here. To estimate 
the stratification profile of the iron abundance we fit the 
data with a straight line using a least-square algorithm. We 
have also checked the dependence of each element's abun- 
dance with respect to the lower level excitation potential Ei 
for the sample of analyzed lines (see as an example Fig. [31). 



Slopes of the linearly approximated dependence of each el- 
ement's abundance with respect to log-E^ and logrsooo are 
given in Table [4] for the case of zero microturbulent veloc- 
ity. Table S] shows a clear correlation between the slope of 
the abundance versus logrsooo and log Ei for the elements 
considered in each star. This correlation strongly suggests 
that the obtained stratification is physical and not due to 
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b) 




-3 

logTsooo 

Fig. 6. The same as in Fig. [TJ but for Fell lines extracted for (a) B84 and (b) WF4-3485. 



uncertainties in the description of the temperature profile 
of the stellar atmosphere. 

The detected stratification of iron, obtained in some 
of the stars studied here, is far too large (> 2 dex) to 
be interpreted as the result of measurement errors. Our 
stellar atmosphere models were calculated taking into ac- 
count the enhanced iron abundance and helium depletion 
(Behr . 2003b 1 and the error in estimating the iron abun- 
dances, for a given model atmosphere, is expected to be 
less than 0.2 dex. Meanwhile, Khan & Shulyak (pUg7|) 
have shown that using atmospheric models with varying 
Fe abundance (from one to ten times solar) can modify the 
abundances inferred from observed line profiles by up to 
±0.25 dex. 

From the results in Table U] we can conclude that no 
detectable iron stratification exists in WF4-3485 and B84. 
Two stars (B267 and B279) show quite large slopes of iron 
abundance with respect to both log Ei and log tsooo • The 
corresponding slopes for the stars WF2-2541 and WF2-3085 
are smaller but still statistically significant. The least signif- 
icant of these slopes is that of the abundance versus log Ei 
for WF2-3085 which has only approximately a 3cr value. 

Another factor that must be evaluated that could mimic 
stratification is an error in the effective temperature of the 
underlying atmospheric model used. To verify the potential 
importance of this factor, we calculated the Fe stratification 
with models with an effective temperature higher and lower 
by 1000 K from the T^g values listed in Table H We also 
calculated the stratification with models assuming solar He 
and Fe abundances to evaluate the effect of this change 
on the inferred stratification profiles. Our simulations show 
that models with higher and models with solar abun- 
dances usually slightly decrease the slope of the abundance 
with respect to log Ei, while models with lower Tcs increase 
this slope for all elements. In the case of iron in WF4-3085, 
the model with higher T^s and the model with solar abun- 
dances result in a negligible slope of its abundance versus 
log Ei. Therefore, the results are not sufficient to report 
unambiguous detection of vertical stratification of iron in 
WF4-3085. The slopes for iron found for the other three 
stars (B267, B279 and WF-2541) are stiU significant in aU 
these models, and stratification is therefore confirmed. 



The two BHB stars (B84 in M15 and WF4-3485 in M13) 
which do not show any signs of stratification of their iron 
abundances are not different from the other stars with re- 
gards to their stellar atmosphere characteristics and rota- 
tion (see Table [5]). However, these two stars possess very 
small titanium and phosphorus abundances as compared 
to the other stars studied here. Their spectra have very 
weak or invisible Ti ii and P il lines. Their iron abundance 
is also close to the solar value, at least for ^ = Okms~^ (see 
Table O. 

We observe an upturn in the iron abundance obtained 
assuming ^ = Okms""'^ at low optical depths for some of the 
stars studied (e.g. Figs. [1^, [5^, S^, [5^). The points on the 
figures corresponding to the upturn feature are the strong 
Fell lines and are not taken into account during the least- 
square fit of the data. To estimate the influence of micro- 
turbulence we have performed an abundance analysis with 
^ = 2kms^^, which leads to elimination of the upturn (see, 
for example, Fig.[T^). Performing a set of simulations with 
different microturbulent velocities, we have determined the 
minimum value of the microturbulence (^mm) for which 
the upturn disappears for each star studied. These values 
are reported in Table [5] As the inclusion of microturbulent 
velocity amplifies the iron stratification, the stratification 
obtained for iron in the first four stars studied with ^ = 
Okms~^ can be considered to be a lower limit to the pos- 
sible stratification in these stars. A more detailed study 
including abundance stratification in the calculation of the 
synthesized spectra might shed some light on this strange 
behaviour, but this is outside of the scope of the present 
paper. 

Figure [T}d indicates that B267 shows a possible signa- 
ture of vertical stratification of the titanium abundance. 
The respective slopes of Ti abundance versus log Ei and 
log T5000 are significantly higher than zero for this star and 
are statistically significant (see Table |4]). For B279, we can 
not confidently conclude that stratification exists since the 
slope of the abundance with respect to log Ei becomes weak 
when using either solar abundances or assuming a T^s in- 
crease of 1000 K. The slope of Ti abundance with respect to 
logTsooo also becomes weak in the model assuming a Toff 
increase of 1000 K. It should be noted that the Ti lines 
in B279 do not sample a large portion of the atmosphere 
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Table 5. Characteristics of our sample of hot BHB stars obtained from spectral simulations. 



Cluster/Star 


vsin I 








[Fe/HJ 






[P/H| 






km 8 


km s 


, _ 1 

km s 


^=0 km s"^ 


5=2 km s ^ 


n 


5=0 km s"^ 


5=2 km s ^ 


n 


M13/WF2-2541 


2.4±0.5 


-257.3±0.5 


0.0+0.3 


-0.12+0.20 


-0.19+0.26 


53 


+1.51+0.11 


+ 1.46+0.10 


16 


M13/WF4-3085 


3.4±0.8 


-255.4±0.5 


0.5±0.3 


+0.72+0.29 


+0.59+0.37 


117 


+1.84+0.13 


+ 1.84+0.15 


12 


M13/WF4-3485 


3.0±0.9 


-246.9±0.9 


1.5±0.3 


-0.04+0.50 


-0.17+0.27 


50 








M15/B84 


4.0±1.3 


-108.3±0.4 


0.0+0.3 


-0.14+0.10 


-0.35+0.24 


25 






















[Ti/H] 




M15/B267 


6.2±1.2 


-114.7±0.7 


0.5+0.3 


+0.45+0.37 


+0.21+0.51 


46 


+0.09+0.20 


-0.08+0.24 


16 


M15/B279 


4.7±1.6 


-104.5±0.6 


1.0+0.3 


+0.42+0.34 


+0.20+0.40 


71 


-0.44+0.10 


-0.54+0.11 


11 



and this fact renders the detection of stratification more 
difficult. 

It is clear that phosphorous shows no clear signs of strat- 
ification in WF2-2541. For WF4-3085 we cannot make firm 
conclusions concerning phosphorous stratification since its 
abundance variation in the range of optical depths under 
consideration is only approximately 0.4 dex. Also, the slope 
of its abundance with respect to log Ei is not statistically 
significant. 

5. Mean abundances and velocities 

The mean photospheric abundances derived for the selected 
BHB stars from an analysis of the available iron, titanium 
and phosphorus lines are reported in Table [5l For each 
chemical element, three columns are shown containing the 
mean abundance for simulations with microturbulent ve- 
locity ^=0 and 2 km s~^ respectively, and the number of 
analyzed line profiles. The reported uncertainties are equal 
to the standard deviation calculated from the results of in- 
dividual line simulations for all lines considered. 

The derived heliocentric radial and projected rotation 
velocities as well as microturbulent velocities (see Table [5]) 
are generally in agreement with the data provided by Behr 
(j2003b,) for these stars. Some inconsistencies are found 
between wsini values obtained for WF2-2541 (wsini — 
O.Otol^^ km s-i) and B279 {vsini = 5.92ti;^g km s"!), 
but they are within the respective error bars given by Behr 
(|2003b|) . The precision of our velocity estimates appears to 
be comparatively higher because of a more stringent selec- 
tion of analysed line profiles. Our measurements of the aver- 
age abundance of iron, titanium and phosphorus also agree 
with the corresponding values published by Behr (,2003bj ), 
taking into account the error bars and the value of the ap- 
plied microturbulent velocity. Small differences between our 
results and abundances from the aforementioned paper are 
not surprising. Our average abundance is calculated with 
the individual abundances obtained from each line, while 
Behr (|2003bp fitted the whole spectrum of each element to 
obtain its abundance. 

6. Discussion 

In this paper we continue our attempts to detect vertical 
abundance stratification in the atmospheres of BHB stars. 
After the report by Bonifacio et al. (|l995p of the vertical 
stratification of helium in the atmosphere of Feige 86, it be- 
came clear that the abundances of other chemical species 
may also be stratified. We devoted special interest to iron 
because hot BHB stars usually have an enhanced iron abun- 
dance (e.g. Behr l2003bp . suggesting that this element may 



be strongly affected by diffusion. Analysing the spectra of 
another hot BHB star HD 135485 (Khalack et al. [2007|) we 
did not find direct evidence of iron stratification, but re- 
vealed strong signatures of sulfur depletion in the deeper 
atmospheric layers. However, HD 135485 is different from 
the other BHB stars in that its spectrum shows evidence 
of helium enrichment (in comparison with the solar abun- 
dance), while in the atmospheres of the other BHB stars 
helium is depleted. 

Therefore, we directed our attention to BHB stars where 
the iron abundance is near the solar abundance or en- 
hanced, and helium is depleted. The results obtained ar- 
gue that at least three stars (B267 and B279 in M15 and 
WF2-2541 in M13) show clear signatures of vertical strat- 
ification of their iron abundance, while for WF4-3085 the 
results are suggestive, but not conclusive. The other two 
stars studied here (B84 in M15 and WF4-3485 in M13) do 
not show stratification of iron and their averaged iron abun- 
dance is close to solar (but is enhanced in comparison with 
its cluster value). B267 shows also a signature of vertical 
stratification of titanium (see Fig. [T|d). 

Since our simulations show that the turnup feature ob- 
served in the iron stratification profile is strongly dependent 
on microturbulent velocity, the value of the abundance at 
low optical depths is uncertain. Of course, if the theoretical 
framework supposes that these abundance gradients are due 
to atomic diffusion, microturbulence should be weak since a 
stable atmosphere is needed for diffusion to be dominant. It 
should be noted that for corresponding optical depths the 
abundance profile of iron is similar in the three BHB stars 
that exhibit stratification. The reason that the other two 
stars in our study do not show clear signs of iron stratifica- 
tion (B84 and WF4-3485) might be related to evolutionary 
effects or the presence of other competing hydrodynamical 
processes. The absence of Ti ii and P ii lines in their spectra 
might be evidence of this. 

In conclusion, the results shown here add to the mount- 
ing evidence of the existence of vertical abundance strati- 
fication, and hence atomic diffusion, in the atmospheres of 
BHB stars. 
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